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PREFACE 


This  report  is  published  to  provide  coastal  engineers  with  the  documenta¬ 
tion  of  a  package  of  computer  programs  for  editing,  analyzing,  and  displaying 
beach  profile  survey  data.  This  package,  named  the  Beach  Profile  Analysis 
System  (BPAS),  was  needed  for  the  analysis  of  a  large  data  bank  of  field  and 
laboratory  profile  surveys.  The  work  was  carried  out  under  the  U.S.  Army 
Coastal  Engineering  Research  Center's  (CERC)  Beach  Profile  Studies  work  unit. 
Shore  Protection  and  Restoration  Program,  Coastal  Engineering  Area  of  Civil 
Works  Research  and  Development. 

This  report  (Vol.  VIII),  the  last  of  eight  volumes,  contains  supporting 
appendixes  that  provide  information  on  data  collection  formatting  and  restric¬ 
tions,  computational  procedures,  assumptions,  and  error  messages  not  repeated 
in  the  other  volumes. 

The  report  was  prepared  by  Marilyn  V.  Fleming,  Systems  Analyst,  under  the 
supervision  of  P.  Pierce,  Chief,  ADP  Office,  with  the  assistance  of  Allan  E. 
DeWall,  Geologist,  under  the  supervision  of  C.J.  Galvin,  former  Chief,  Coastal 
Processes  Branch,  and  Mr.  R.P.  Savage,  Chief,  Research  Division. 

Instrumental  insight  concerning  a  previous  version  of  the  Beach  Profile 
Analysis  System  was  provided  by  B.  Sims.  Programing  was  accomplished  by 
M.  Fleming  and  T.  Lawler  with  the  assistance  of  D.  French,  J.  Alquist, 
R.  Hylton,  and  F.  Wilson. 

The  authors  acknowledge  the  helpful  discussions  and  review  comments  of 
Drs.  C.  Everts,  C.  Galvin,  R.  Hallermeier,  and  C.  Vincent,  and  W.  Birkemeier, 
M.  Hemsley,  T.  Lawler,  H.C.  Miller,  B.  Sims,  and  P.  Vitale. 

Technical  Director  of  CERC  was  Dr.  Robert  W.  Whalin,  P.E.,  upon  publica¬ 
tion  of  this  report. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th 
Congress,  approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th 
Congress,  approved  7  November  1963. 


'TED  E.  BISHOP j/ 
Colonel,  Corps  of  Engineers 


Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 

U.S.  customary  units  of  measurement  used  In  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

Inches 

25.4 

millimeters 

2.54 

centimeters 

square  Inches 

6.452 

square  centimeters 

cubic  Inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius 

(C) 

temperature  readings 

from  Fahrenheit  (F)  readings, 

use  formula:  C  ■  (5/9) 

(F  -32). 

To  obtain  Kelvin 

(K) 

readings,  use  formula: 

K  =  (5/9)  (F  -32)  +  273.15. 
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APPENDIX  A 


GLOSSARY  OP  TERMS 

BENCH  MARK  -  A  permanently  fixed  point  of  known  position  and  elevation. 

CONTOUR  -  A  line  of  constant  elevation  along  the  beach  surface. 

CONTOUR  INTERCEPT  -  The  point  defined  by  the  seawardmost  intersection  of  a 
contour  with  a  beach  profile.  On  some  profiles,  there  may  be  more 
than  one  intercept  of  a  given  contour,  (e.g.,  MULTIPLE  INTERCEPT). 

CONTOUR  POSITION  -  The  intersection  of  a  horizontal  plane  and  the  beach 
surface . 

DATUM  INTERCEPT  -  The  seawardmost  point  where  a  profile  crosses  the  VERTICAL 
DATUM  plane  (e.g.,  the  SHORELINE  POSITION). 

DATUM,  VERTICAL  -  The  zero  elevation. 

INPUT  -  The  vertical  datum  to  which  survey  input  data  are  referenced. 

OUTPUT  -  The  vertical  datum  to  which  the  data  represented  in  output  plots 
and  tables  are  referenced. 

DATUM,  HORIZONTAL  -  The  zero  distance. 

INPUT  -  The  horizontal  datum,  usually  BENCH  MARK,  to  which  survey  input 
data  are  referenced. 

OUTPUT  -  The  horizontal  datum  to  which  the  data  represented  in  output  plots 
and  tables  are  referenced  (e.g.,  SHORELINE  POSITION  during  REFERENCE 
SURVEY,  dune  crest  position,  etc.). 

DEFAULT  VALUE  -  The  value  automatically  assigned  to  a  variable  when  no  value  is 
assigned  by  the  user. 

DISTANCE  -  The  horizontal  coordinate  of  a  point  on  a  beach  PROFILE,  measured 
positively  seaward  from  a  fixed  point  on  the  PROFILE  LINE  (e.g.,  BENCH 
MARK). 

ELEVATION  -  Vertical  coordinate  of  a  point  on  a  beach  PROFILE  measured 
positively  upward  from  a  known  VERTICAL  DATUM. 

EXTRAPOLATED  DATUM  -  The  SHORELINE  position  when  it  has  been  extrapolated  using 
the  last  two  surveyed  points. 

INPUT  UNITS  -  Units  of  measurement  in  which  the  survey  input  data  are  recorded. 

LANDWARD  BOUNDARY  -  The  distance  identifying  the  landward  end  of  the  PROFILE. 

LOCALITY  -  A  unique  place  where  survey  data  are  located  (e.g.,  Atlantic  City, 
Bodie  Island,  etc.). 


LOCALITY  CODE  -  Two-character  representation  of  the  locality. 


LOCALITY  DESCRIPTION  -  User-supplied  name  of  the  locality  or  a  description 
of  the  data. 

MULTIPLE  INTERCEPT  -  A  point  where  a  specific  CONTOUR  crosses  the  PROFILE  that 
is  in  addition  to  the  seaward  intercept. 

OUTPUT  UNITS  -  Units  of  measurement  of  data  displayed  on  output  plots  and 
tables. 

PROFILE  -  The  curve  defined  at  the  intersection  of  the  ground  surface  with  a 
vertical  plane  inserted  perpendicular  to  it. 

PROFILE  AREA  -  The  area  bounded  above  by  a  PROFILE,  landward  by  a  vertical 
line,  below  by  a  horizontal  line  and  seaward  by  a  vertical  line. 

PROFILE  ENVELOPE  -  The  area  bounded  by  the  maximum  and  minimum  elevations 
found  at  distances  along  a  PROFILE  LINE. 

PROFILE  LINE  -  The  line  defined  by  two  fixed  points  or  by  one  fixed  point 
and  a  known  angle,  along  which  surveyors  measure  distances  and 
elevations  to  define  a  beach  PROFILE. 


PROFILE  LINE  NUMBER  -  The  unique  number  assigned  to  each  PROFILE  LINE. 

REFERENCE  SURVEY  -  The  survey  of  a  PROFILE  LINE  which  defines  a  standard  value 
to  which  subsequent  values  are  referenced  (e.g.,  SHORELINE  position  or  UNIT 

VOLUME) . 


SEAWARD  BOUNDARY  -  The  distance  identifying  the  seaward  end  of  the  PROFILE. 
UNIT  VOLUME  i3  computed  from  this  point  landward. 

SHORELINE  -  The  seawardmost  intercept  of  the  VERTICAL  DATUM  plane  with  the 
beach. 


SPATIAL  MEAN  -  The  mathematical  average  of  a  set  of  variables  collected  from  a 
number  of  PROFILE  LINES  during  a  single  SURVEY. 

SURVEY  -  An  event  during  which  one  or  all  of  the  PROFILE  LINES  at  a  LOCALITY 
are  surveyed. 

SURVEY  NUMBER  -  The  unique  number  assigned  to  each  set  of  profile  data  surveyed 
during  the  same  time. 


TEMPORAL  MEAN  -  The  mathematical  average  of  a  set  of  variables  collected  during 
a  number  of  SURVEYS  at  a  single  PROFILE  LINE. 


UNIT  VOLUME  -  The  product  of  a  cross-sectional  area  and  a  unit  of  length  per¬ 
pendicular  to  the  area.  Given  as  units  of  volume  per  length  of  shoreline 
(e.g.,  m’/m  or  yd^/ft). 


UNIT  VOLUME  BY  CONTOUR  -  Unit  volume  computed  within  slices  bounded  by 
specific  CONTOURS. 
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UNIT  VOLUME,  ZERO  -  The  unit  volume  which  establishes  the  standard  for 
computing  change  in  unit  volume  during  SURVEYS  of  a  PROFILE  LINE. 
This  volume  is  subtracted  from  volumes  found  during  other  SURVEYS. 

UNIT  VOLUME,  CHANGE  -  The  difference  between  unit  volumes  measured  during 
two  SURVEYS. 
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SYMBOLS  AND  DEFINITIONS 


A 

C 

c 

d 

H  -j.  • 
i 

j 

k 

M^» 

M— t 

m 

n 

P 

Q 

r 

S 

s 

Tt;. 

t 

u 

V 

w 

XLAND 


an  area 

conversion  factor  to  convert  area  to  unit  volume 
elevation  of  a  contour 

distance  from  a  profile  line  to  the  adjacent  one 
hours  elapsed  since  the  time  of  first  survey 

integer  identifying  position  (i=l  for  landwardmost ,  of  coordinate- 

pair  in  the  array  defining  a  surveyed  profile  line 

number  of  surveys  of  a  locality 

number  of  surveys  for  which  a  value  can  be  computed  or  defined 

integer  identifying  each  survey 

months  elapsed  since  time  of  first  survey 

months  elapsed  from  the  beginning  of  the  year  of  a  particular  survey 
integer  identifying  each  profile  line  at  a  locality 
number  of  profile  lines  at  a  locality 
a  selected  distance  on  a  profile  line 

the  number  of  distance  and  elevation  pairs  defining  a  survey  of  a 
profile  line 

correlation  coefficient 

slope  of  a  line 

subscript  identifying  a  mean  as  spatial 

zero  time,  the  time  the  earliest  survey  in  the  survey  data  file  was 
taken 

subscript  identifying  a  mean  as  temporal 

distance  along  the  beach  from  the  end  profile  lines  for  which  the 
computed  distance  to  a  contour  may  be  considered  valid 

unit  volume  of  a  profile 

weighting  factor  used  to  compute  spatial  means 
landward  distance  boundary  for  area  computations 
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SYMBOLS  AND  DEFINITIONS— Continued 
XSEA  seaward  distance  boundary  for  area  computations 
X,  Y  independent  and  dependent  variables  used  in  regression  analysis 
x  surveyed  or  computed  distance  coordinate  on  a  profile 

YMAX  upper  elevation  boundary  for  area  computations 

YMIN  lower  elevation  boundary  for  area  computations 

y  surveyed  or  computed  elevation  coordinate 

3  intercept  of  regression  line 

3  regression  coefficient 

a  standard  deviation 
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APPENDIX  B 


COMPUTATIONAL  PROCEDURES 

I.  EXTRAPOLATION  OF  SHORELINE  POSITION 

Since  some  records  may  not  be  included  in  certain  computations  because  the 
profiles  they  describe  do  not  extend  to  the  shoreline,  the  user  may  decide  that 
the  distance  to  the  shoreline  should  be  extrapolated  and  added  a3  a  surveyed 
point.  When  allowing  extrapolation,  the  user  is  assuming  that  the  beach  slopes 
evenly  from  a  selected  minimum  elevation  to  the  shoreline.  The  criteria  for 
extrapolation  are  as  follows: 

(a)  The  elevation  of  the  last  surveyed  point  must  not  exceed  the 
selected  minimum  elevation. 

(b)  The  slope  of  the  line  connecting  the  last  two  surveyed  points 
must  be  negative  (beach  is  sloping  seaward). 

Using  the  last  two  surveyed  points  in  the  record,  (*„•*  >'/ )  and  (xz+r  W 
and  the  elevation  at  the  shoreline,  0,  the  equation  of  the  line  connecting  the 
last  two  points  is  solved  for  the  distance  to  the  0  elevation: 


>-i+l 


>7+2  (xz  -  xi+l} 

{n  - 


(B-1 ) 


This  point  is  added  to  the  record  as  the  last  surveyed  point  and  a  flag  set  for 
the  record  indicating  that  the  shoreline  was  extrapolated.  All  output  results 
computed  using  extrapolated  data  are  flagged.  For  all  further  computations, 
the  record  is  treated  as  if  it  had  a  shoreline. 

II.  COMPUTING  DISTANCES  TO  GIVEN  ELEVATIONS 

This  computation  is  utilized  in  the  BPAS  routines  to  interpolate  the 
distance  to  the  shoreline  and  in  any  instance  when  it  is  necessary  to  find  the 
distance  at  a  specific  elevation,  but  its  most  extensive  application  is  in  the 
analysis  module  SURVY2  which  computes  and  displays  distance  and  changes  in 
distance  to  selected  contours.  To  find  the  seawardmost  distance  to  these 
contours,  each  surveyed  point,  beginning  with  the  seawardmost  and  proceedings 
landward,  is  tested  until  the  elevation  of  the  desired  contour  is  either  the 
same  as  the  elevation  of  a  surveyed  point  or  falls  between  the  elevations  of 
two  adjacent  surveyed  points.  When  the  elevation  is  the  same,  the  distance  is 
as  surveyed  and  no  further  computation  is  required.  Otherwise,  the  equation  of 
the  line  through  the  two  adjacent  points,  (*.,  y^) ,  (x^+1>  >7+1),  is  solved  for 
the  distance  coordinate  at  the  desired  elevation,  c: 


(c 

XC  =  H+1  +  — 


Vi+1 ) 

"(yi  -  yi+1) 


(B-2) 
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For  each  contour,  a  flag  is  set  to  indicate  whether  the  distance  was  surveyed 
or  interpolated,  whether  the  computation  required  the  use  of  an  extrapolated 
shoreline  position,  or  whether  the  desired  contour  was  undefined  by  the  survey 
data. 


At  times,  a  contour  may  be  defined  on  a  profile  more  than  once,  such  as 
when  a  dune  or  a  sandbar  is  surveyed  (Fig.  B-l).  Provision  has  been  made  to 
compute  and  display  the  distance  to  up  to  10  of  these  multiple  intercepts  if 
such  is  desired.  Unless  otherwise  indicated,  however,  only  the  distance  to  the 
seawardmost  will  be  computed.  In  further  computations  involving  contour 
positions,  the  seawardmost  distance  is  used. 


Figure  B-l.  Multiple  contour  intercepts.  The  contours  C3 ,  C4,  and  C5 
all  have  "multiple  intercepts"  on  the  defined  profile. 

III.  DETERMINING  THE  SLOPE  OF  A  PROFILE  AT  THE  SHORELINE 

In  the  BPAS  routines,  the  slope  at  the  shoreline  of  the  profile  defined  by 
each  record  i3  computed.  When  the  shoreline  position  is  interpolated,  the 
slope  of  the  line  through  the  two  points  between  which  it  falls  is  used.  If 
the  shoreline  is  defined  by  a  surveyed  or  extrapolated  point,  the  slope  of  the 
line  through  the  shoreline  and  the  next  landward  point  is  used: 


S 


(B-3) 


IV.  COMPUTING  MEAN  DISTANCE— TEMPORAL  ' 

The  temporal  mean  is  a  simple  mean  used  to  determine  the  mean  position  of  a 
contour  at  a  profile  line  considering  all  surveys  taken  at  that  line.  It  is 
used  to  compute  the  mean  shoreline  position  in  the  BPAS  routines,  when  such  is 
selected  as  output  horizontal  datum,  and  in  the  BEACH  analysis  module.  In 
analysis  module  SURVY2,  the  mean  position  of  selected  contours  at  each  profile 
line  is  computed  using  this  method.  The  mean  is  determined  by  adding  the 
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distance  to  the  contour  for  each  survey  of  the  profile  line  and  dividing  by  the 
number  of  surveys  during  which  the  distance  could  be  computed.  Surveys  during 
which  the  distance  could  not  be  defined,  because  the  contour  did  not  exist,  are 
ignored . 

Letting 

<  be  a  consecutive  integer  identifying  each  survey  of  the  profile  line, 

xs  be  the  distance  to  the  contour  c  during  survey  , 

V 

x  >  be  0  for  surveys  during  which  the  position  of  contour  c  cannot  be 
«.  defined, 

j  be  the  total  number  of  times  the  profile  line  is  surveyed,  and 

k  be  the  number  of  surveys  during  which  an  x,  can  be  computed. 

The  mean  contour  position  for  contour  c  at  each  profile  line: 


(x.). 


(B-4) 


V.  COMPUTING  MEAN  DISTANCE-SPATIAL 

Analysis  module  SURVY2  is  the  only  module  which  utilizes  this  procedure. 

It  yields  the  mean  position  of  the  selected  contours  for  a  given  survey.  The 
position  of  each  contour  on  all  profile  lines  surveyed  is  considered  and  a 
weighted  mean  is  computed.  The  weighting  factor  for  each  profile  line  is 
unique  and  based  on  the  following  assumptions: 

(a)  The  profile  lines  are  parallel. 

(b)  Each  profile  is  representative  of  the  beach  for  a  distance 
halfway  to  the  adjacent  ones;  the  user  must  provide  an  appropriate 
representative  distance  for  the  end  profile  lines. 

The  U3er  provides  the  distance  between  the  adjacent  profile  lines;  if  none  are 
supplied,  the  lines  are  assumed  to  be  equidistant  and  a  simple  mean  is 
computed.  Using  the  profile  spacing  provided,  the  weighting  factor  is  computed 
for  each  profile  line  by  adding  half  the  distance  to  the  profile  line  to  its 
left  to  half  the  distance  to  the  one  on  its  right  (Fig.  B-2).  The  distance  to 
each  contour  at  each  profile  line  is  multiplied  by  the  weighting  factor  for 
that  line  and  the  results  added  for  all  the  lines.  This  sum  is  divided  by  the 
sum  of  all  the  weighting  factors  to  determine  the  mean  spatial  position  for 
each  contour. 
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Profile  Lme  I 


Profile  Line  2 


Profile  Line  3 


Profile  Line  4 


Profile  Lme  5 


WATER 


Profile  Line  n-l 


Profile  Line  n 


LAND 


WATER 


Figure  B-2.  Distances  used  to  determine  the  weighting  factor. 


Letting 

m  be  a  consecutive  integer  identifying  each  profile  line  at  the 

locality, 

n  be  the  number  of  profile  lines  at  the  locality, 

U|,  \in  be  the  representative  distances  for  the  end  profile  lines,  and 

dm  be  the  distance  from  profile  line  m  to  m  +  1  (Fig.  B-2). 

The  weighting  factors  for  the  end  profile  lines 


wi  =  ui  +  T 


(B-5) 


,  u-1 

w  =  u  +  — ^ — 

n  >:  2 


(B-6) 
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The  weighting  factor  for  each  internal  profile  line 


"m 


+ 


(B-7) 


and  letting 


w  m 


be  the  distance  to  the  seawardmost  intercept  of 

be  0  for  profile  lines  not  surveyed  or  at  which 
contour  c  is  not  defined. 


The  mean  position  of  the  contour  c  during  each  survey 


contour  c,  and 
the  position  of 


(*e)s 


l 

m=l 


w  x 

m  a 


n 


£  w K 

m=l 


(B-8) 


It  should  be  noted  that  when  a  contour  cannot  be  defined  at  a  profile  line,  it 
is  as  if  the  profile  line  and  the  part  of  the  beach  represented  by  it  do  not 

exist  (Fig.  B-3). 


WATER 


Figure  B-3. 


WATER 


Spatial  mean  contour  position  when  the  contour  fur  which 
the  spatial  mean  is  being  computed  is  not  defined  at 
profile  line  2. 
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VI.  COMPUTING  MEAN  DISTANCE— SPATIAL-TEMPORAL 

Used  exclusively  by  analysis  module  SURVY2,  the  spatial-temporal  mean  is 
computed  to  give  an  average  profile  at  a  locality  as  defined  by  all  surveys 
taken  at  that  locality.  It  is  the  temporal  mean  of  the  position  of  each 
selected  contour  computed  using  the  spatial  mean  of  the  contour  found  during 
each  survey. 

Letting 

£  be  a  consecutive  integer  identifying  each  survey  at  a  locality, 

0x^)3 i  be  the  mean  position  of  the  contour  c  at  the  locality  during 

survey  f , 

be  0  for  surveys  during  which  it  cannot  be  computed, 
j  be  the  number  of  surveys  at  the  locality,  and 

k  be  the  number  of  surveys  at  the  locality  during  which  (xe)s?  was 

defined. 

The  mean  contour  position  for  contour,  c 


j 

£*i 

k 


(B-9) 


VII.  COMPUTING  ELEVATIONS  AT  FIXED  DISTANCES 

The  procedure  used  for  this  computation  is  quite  similar  to  the  one  used  to 
compute  distances  at  specific  elevations.  The  major  difference  is  that  there 

is  usually  only  one  elevation  at  a  specific  distance  during  a  given  survey  of  a 
profile  line.  When  there  is  a  scarp  defined  by  two  elevations  at  a  distance, 
the  elevation  with  the  higher  subscript  is  used.  Although  used  whenever  it  is 
necessary  that  the  elevation  at  a  specific  distance  be  known,  this  procedure 
is  mainly  used  in  the  analysis  module  ELVDIS  to  compute  and  display  elevations 
and  elevation  changes  and  the  maximum  and  minimum  elevations  at  fixed  distances 
along  each  profile  line. 


Beginning  at  the  seawardmost  point,  each  coordinate  pair  is  tested  until 
the  fixed  distance  under  consideration  is  either  equal  to  the  distance  coordi¬ 
nate  of  a  surveyed  point  or  found  to  lie  between  two  adjacent  surveyed  points. 
Of  course,  when  the  fixed  distance  is  at  a  surveyed  point,  the  corresponding 
elevation  is  known.  If  not,  the  equation  of  the  line  through  the  two  adjacent 
points  (x  . ,  y. ),  (x.  ,  x.  ),  is  solved  for  the  elevation  at  the  fixed 

...  u  I*  Z-*r  1  2-  '  1 

distance,  p: 
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(B-10 ) 


yp  =  n+i 


(p  -  h+i]^l  -  n+i] 

( h  -  H+0 


For  each  fixed  distance  a  flag  is  set  to  indicate  whether  the  associated  eleva¬ 
tion  was  surveyed  or  found  through  normal  interpolation,  whether  the  computa¬ 
tion  required  the  use  of  an  extrapolated  shoreline,  or  whether  the  fixed  dis¬ 
tance  was  undefined  by  the  survey  data. 

VIII.  DETERMINING  THE  MAXIMUM  AND  MINIMUM  ELEVATIONS 
AT  FIXED  DISTANCES  ALONG  A  PROFILE  LINE 

Using  the  procedure  defined  above,  elevations  are  found  at  fixed  distances 
along  the  profile  line,  compared,  and  the  maximums  and  minimums  saved.  The 
profile  envelope  thus  formed  is  composed  of  the  maximums  and  minimums  for  all 
surveys  at  the  profile  line  or,  if  desired,  only  for  surveys  during  each  year 
covered  by  the  data.  In  the  latter  case,  there  will  be  an  envelope  for  every 
year  during  which  there  was  at  least  one  survey  at  the  profile  line. 

IX.  COMPUTING  UNIT  VOLUME 

The  analysis  modules  in  the  BP AS  using  unit  volume  are  BEACH  and  VOLCTR. 

The  unit  volume  is  actually  a  volumetric  representation  of  an  area  determined 
by  computing  the  area  within  boundaries  defined  by  an  upper  and  lower  eleva¬ 
tion,  a  landward  and  seaward  distance,  and  the  profile.  This  area  is  extended 
laterally  along  one  unit  length  of  the  coast  to  give  a  volumetric  representa¬ 
tion  of  the  beach  at  the  profile  line  (Fig.  B-M).  If  the  lower  elevation 
boundary  is  less  than  zero,  two  unit  volumes  are  computed:  one  for  the  area 
below  the  zero  elevation  and  one  for  the  area  above. 

All  unit  volume  computations  are  performed  in  a  single  subroutine.  The 
coordinate  pairs  describing  the  profile,  the  upper  and  lower  elevation  bound¬ 
aries,  and  the  landward  and  seaward  distance  boundaries  are  passed  to  this  sub¬ 
routine.  The  subroutine  computes  the  area  within  the  given  boundaries,  two 
areas  when  the  area  below  the  zero  elevation  is  required,  converts  area  to  unit 
volume,  and  returns  the  unit  volume  to  the  calling  routine.  The  necessary 
boundaries  for  the  area  computation  are  determined  by  the  calling  routine. 

1 .  Determining  Boundaries 

When  discussing  boundaries,  reference  is  made  to  the  previous  description 
in  Volume  I  (Sec.  IV)  of  the  three  types  of  profiles  which  can  be  processed 
by  the  BPAS. 

a.  Analysis  Module  VOLCTR — Total  Unit  Volume  and  Unit  Volume  by  Contour. 

In  analysis  module  VOLCTR,  two  types  of  unit  volume  are  computed:  total  unit 
volume  and  unit  volume  by  contour.  The  total  unit  volume  in  VOLCTR  is  computed 
with  as  much  of  the  survey  data  as  possible  for  each  survey  of  a  profile  line 
(Fig.  B-5)-  There  are  no  user-supplied  or  computed  upper  and  lower  elevation 
boundaries  or  landward  and  seaward  distance  boundaries.  For  type  1  and  f  pro¬ 
files,  the  part  of  the  area  above  the  zero  elevation  within  the  following 
boundaries  is  computed. 
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Zero  Elevotion 


Profile  Type  2 


Profile  Type  3 

Figure  B-5.  Boundaries  for  total  unit  volume  computations — VOLCTP . 
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Landward — landwardmost  distance 

Seaward - seawardmost  intercept  of  the  zero  elevation 

Upper - profile  surface 

Lower - zero  elevation 


For  type  3  profiles,  the  landward  and  upper  boundaries  remain  unchanged.  The 
others  are: 


Lower - elevation  of  last,  not  lowest  surveyed  point 

Seaward - distance  to  last  surveyed  point 


The  purpose  of  the  unit  volume  by  contour  computations  in  analysis  module  VOLCTR 
is  to  compare  unit  volume  changes  within  specific  contour  intervals  from  one 
survey  of  a  profile  line  to  the  next;  the  overall  boundaries  used  are  therefore 
determined  based  on  common  area  between  the  two  surveys.  Although  these  com¬ 
parisons  are  made  only  for  the  part  of  the  profiles  above  the  zero  elevation, 
it  is  not  necessary  that  a  profile  extend  to  the  zero  elevation  in  order  to 
qualify  for  these  unit  volume  computations.  The  comparison  is  made  for  all 
contour  intervals  which  the  consecutive  surveys  of  the  profile  line  have  in 
common  (Fig.  B-6).  The  horizontal  segments  within  which  unit  volume  changes 
are  compared  are  consistent.  The  contour  defining  the  bottom  of  the  lowest 
segment  is  zero  or  some  user-selected  minimum  greater  than  zero  and  each  segment 
is  the  same  width,  normally  one  unit.  The  overall  boundaries  are  computed  and 
then  the  unit  volume  within  the  appropriate  segments  is  computed  and  the  change 
within  each  segment  determined  by  subtracting  the  unit  volume  found  during  the 
previous  survey  from  the  one  for  the  current  survey. 


The  upper  and  lower  boundaries  for  the  volume  computations  do  not  always 
fall  directly  on  a  contour  defining  the  limit  of  one  of  the  horizontal  segments. 
When  this  happens  on  the  upper  part,  there  is  no  problem  since  the  profile 
surface  serves  as  the  upper  boundary  for  the  volume  computation.  However,  when 
it  happens  at  the  lower  boundary,  it  is  necessary  to  compute  the  width  of  a 
partial  segment  for  the  first  comparison.  Unit  volumes  within  these  partial 
intervals  are  output  for  TABLE10  (Vol.  I,  Fig.  37)  and  PL0T11  (Vol.  1,  Fig. 

39) *  and  are  not  output  for  TABLE11  (Vol.  I,  Fig.  38)  or  PL0T12  (Vol.  I,  Fig. 

10). 


The  boundaries  for  the  overall  computation,  based  on  the  data  from  the  two 
surveys  of  the  profile  line  to  be  compared,  are  as  follows: 


Landward — the  landwardmost  distance  in  common  to  both  surveys  of  the 
profile  line 

Seaward - the  seawardmost  intercept  of  the  lower  bound 


Upper - the  profile  surface  or  a  user-defined  upper  contour, 

whichever  is  lower 
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Landward  Bound  |  |  Landword  Bound  \  |  ^  |  |  Landward  Bound 


(1) 


! 


Determine  common  boundaries 


)  Divide  profiles  into  horizontal  segments 


Survey  £ 


(4)  Change  in  unit  volume  at  segment  1  from  survey  ■  to  survey  ■  +1  at  the 
profile  line  is  unit  volume  for  segment  1  during  survey  •  -*-1  minus  unit  volume 
for  segment  1  during  survey 

Figure  B-b.  Unit  volume  by  contour  (consecutive  surveys 
at  the  same  profile  1 i no ) . 
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Lower 


•the  zero  elevation,  a  user-defined  lower  contour  or  the 
greater  of  the  elevations  at  the  seawardmost  point  of  the 
two  surveys,  whichever  of  the  three  is  highest 


The  boundaries  for  each  segment  which  falls  within  the  common  area  are  passed 
separately  to  the  subroutine  which  computes  and  returns  the  unit  volume  for 
the  segment.  The  calling  routine  computes  the  change. 

b.  Analysis  Module  BEACH — Total  Unit  Volume  Above  and  Below  Vertical 
Datum.  The  unit  volume  computations  in  analysis  module  BEACH  determine  changes 
in  unit  volume  at  the  various  profile  lines.  The  boundaries  for  the  unit 
volume  computations  are  therefore  determined  by  all  surveys  taken  at  specific 
profile  lines.  Some  surveys  of  a  profile  line  may  be  eliminated  from  computa¬ 
tions  when  certain  criteria  are  not  met  and,  since  there  are  more  requirements 
imposed  for  unit  volume  below  vertical  datum  computations,  it  is  possible  that 
for  a  given  survey,  the  unit  volume  above  the  vertical  datum  will  be  computed 
while  the  unit  volume  below  will  not. 

The  criteria  for  unit  volume  above  the  vertical  datum  computations  are  as 

follow. : 


(l)  The  survey  must  describe  a  type  1  or  a  type  C  profile  (Vo! .  I, 
fee.  IV). 


(2)  The  distance  to  the  landwardmost  point  must  be  less  than  or  equal 
to  the  input  horizontal  datum  (zero).  This  criteria  may  be  changed  if  the 
user  wishes  to  supply  a  minimum  beginning  distance.  If  the  reference  unit 
volume  is  user-supplied,  the  user  must  also  supply  a  landward  boundary  for 
these  computations;  the  profile  must  begin  no  farther  seaward  than  this 
boundary. 

The  boundaries  for  the  unit  volume  above  vertical  datum  computations  in 
analysis  module  BEACH  are: 

Landward — if  not  user-supplied,  the  landwardmost  point  in  common 
to  all  surveys  of  the  profile  line  which  meet  the 
necessary  criteria 

Seaward - seawardmost  intercept  of  the  zero  elevation.  It  is 

probable  that  this  will  be  at  a  different  point  for 
each  survey  of  the  profile  line 

Upper - the  profile  surface  or  a  user-selected  elevation, 

whichever  is  lower 

Lower - zero  elevation 

The  criteria  for  the  unit  volume  below  vertical  datum  computation-.': 

(1)  The  user  rust  specify  that  the  unit  volume  below  vertical 
datum  computations  are  to  be  performed. 

(?)  The  profile  must  be  a  type  1  profile  which  meets  the  cri¬ 
teria  for  unit  volume  above  vertical  datum  computations . 
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(3)  The  user  must  supply  as  a  tolerance  a  distance  seaward  of  the  output 
horizontal  datum.  If  the  last  surveyed  point  is  less  than  this  tolerance,  the 
survey  of  the  profile  line  is  eliminated  from  further  processing  for  unit  volume 
below  the  vertical  datum.  When  the  reference  unit  volume  below  is  supplied  ly 
the  user,  tolerance  for  each  profile  line  must  also  be  supplied. 

The  boundaries  for  the  unit  volume  below  vertical  datum  computations  in 
analysis  module  BEACH  are: 

Landward — landward  boundary  used  for  the  unit  volume  above  vertical 
datum  computation 

Seaward - if  not  user-supplied,  the  seawardmost  distance  is  common 

to  all  surveys  of  the  profile  line  which  meet  the  neces¬ 
sary  criteria 

Upper - the  zero  elevation 

Lower - 100  units  or  some  user-supplied  lower  elevation.  If  an 

elevation  is  less  than  the  defined  lower  boundary,  a  message 
is  printed  and  the  lower  boundary  is  changed  to  the  new 
minimum.  This  becomes  the  lower  boundary  for  all  unit 
volume  below  datum  computations  at  the  profile  line  being 
processed  and  for  all  profile  lines  to  be  processed.  It 
does  not  affect  previously  analyzed  profile  lines. 

For  the  unit  volume  below  vertical  datum  computations,  a  vertical  is  dropped 
from  the  profile  elevation  at  the  seaward  boundary  to  the  elevation  of  the 
lower  boundary  (as  shewn  in  Fig.  B— L ) . 

2.  The  Computation. 

As  discussed  earlier,  subroutine  VOLCTR  receives  survey  data  and  boundary 
information,  computes  the  area,  and  converts  it  to  unit  volume.  When  VCLCTF 
receives  the  survey  data,  it  first  examines  the  data  to  determine  if  there  is 
an  area  to  be  computed  below  vertical  datura.  If  there  is,  the  area  for  the 
part  below  is  computed  and  then  the  area  for  the  part  above  is  computed.  When 
the  area  for  the  part  above  is  being  computed,  the  seawardmost  intercept  of  the 
zero  elevation  becomes,  for  computational  purposes,  the  seawardmost  point. 

After  boundaries  are  established,  the  area  is  computed.  Beginning  with 
the  seawardmost  point  and  moving  landward,  each  set  of  adjacent  coordinate 
pairs  is  examined  until  the  seaward  boundary  has  been  crossed.  The  computa¬ 
tions  for  area  begin  here  and  continue  until  the  landward  boundary  is  crossed. 
Each  segment  of  the  profile,  as  defined  by  two  adjacent  coordinate  pairs,  is 
considered  and  the  area  beneath  each  segment  which  falls  within  the  defined 
boundaries  is  computed  and  accumulated.  The  total  area  is  then  multiplied  by 
the  supplied  conversion  factor  to  yield  unit  volume  and  the  result  returned  to 
the  calling  routine.  Figure  B-7  provides  a  schematic  of  the  procedure  and, 
because  the  number  of  possible  paths  prohibits  a  simpler  representation  of 
the  computation,  the  logic  of  the  subroutine  is  provided  as  follows. 
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Fipure  B-7 .  Survey  data  and  boundary  information  received 
by  VOLCTR. 


The  part  of  the  total  area  represented  by  each  line  segment,  from  coordi¬ 
nate  (xj,  )  to  ( x^+_ ,  ) ,  is  called  A^ .  Each  is  broken  down  as  neces¬ 

sary  so  that  it  is  composed  of  rectangles  and  triangles;  up  to  three  components 
may  be  required  to  thus  represent  a  given  A^  (Fig.  F- 8 ) .  These  pieces  will 
be  called  AI ,  A2,  and  A3  and  will  be  reset  to  zero  before  each  is  commuted. 

For  simplicity,  also  let  (XI,  Y1)  represent  (x-,  yz-)  and  (X2,  Y2)  represent 
(x,-+i,  yt‘+j  ) .  These  also  will  be  appropriately  reset  before  the  computation  of 
the  successive  A,-'s.  The  boundaries  are  represented  as  follows. 

YMAX — the  upper  elevation  boundary 

YMIN - the  lower  elevation  boundary 

XSEA — the  seaward  distance  boundary 
XLAND — the  landward  distance  boundary 

Because  of  the  control  in  the  calling  routine,  there  will  never  be  a  case 
where  the  seawardmost  point  is  landward  of  the  seaward  boundary  or  that  the 
landwardmost  point  is  seaward  of  the  landward  boundary.  The  changing  of  the 
values  of  (XI,  Y1)  and  (X2,  Y2)  does  not  affect  the  values  of  (*;,  y/)f 


Beginning  at  the  seawardmost  point,  each  consecutive  set  of  coordinate 
pairs  is  tested  until  all  the  following  conditions  are  met: 

(a)  XI  -  XSEA 

(b)  Y1  YMIN  .  r  Y2  >  YMIN 

(c)  XI  i  X2 

Condition  (a)  indicates  the  area  computations  may  begin  while  conditions  (b) 
and  (c)  assure  that  A  .  /  0. 

Determine  if  any  boundaries  are  crossed  and  make  appropriate  adjustment  if 
so. 

If  XI  <  XLAND 

then  XI  =  XLAND 

Use  equation  ( B— 10)  to  find  the  elevation  at  XLAND. 

Let  this  be  Y1. 

If  X2  •  XSEA 

then  X2  =  XSEA 

Use  equation  (B-10)  to  find  the  elevation  at  XSEA. 

Let  this  be  Y2. 

If  Y1  -  YMIN 

Then  Y1  =  YMIN 

Use  equation  (B-2)  to  find  the  distance  to  YMIN. 

Let  this  be  XI . 


Figure  B-8 .  Component  parts  for  area  computations. 


If  Y2  <  YMIN 

then  Y2  =  YMIN 

Use  equation  (B-2)  to  find  the  distance  to  YMIN. 

Let  this  be  X2. 

If  Y1  <  YMAX  xvJ  Y2  >  YMAX 

then  Y2  =  YMAX 

Use  equation  (B-2)  to  find  the  distance  to  YMAX. 

Let  this  be  X3. 

A1  =  (X2  -  X3)  (Y2  -  YMIN) 

Let  X2  r  X3. 

If  Y1  >  YMAX  Y2  <  YMAX 

then  Y1  =  YMAX 

Use  equation  (B-2)  to  find  the  distance  to  YMAX. 

Let  this  be  X3. 

A1  =  (X3  -  XI)  ( Y2  -  YMIN) 

If  Y1  >  YMAX 

then  Y1  =  YMAX 

If  Y2  >  YMAX 

then  Y2  =  YMAX 

Compute  the  area  under  the  line  segment  and  within  the  pertinent  boundaries 

Y3  =  the  greater  of  Y1  and  Y2 
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Y4  =  the  lesser  of  Y1  and  Y2 
A2  =  (X2  -  XI)  (Y3  -  Y4)/2 
A3  =  (X2  -  XI)  (Y4  -  YMIN) 

A  =  A1  -f  A2  +  A3  (Fig.  B-13) 

The  next  two  coordinate  pairs  are  considered  and  the  next  A;  computed 
until  X2  <  XLAND.  The  required  area  has  then  been  computed. 

Letting 

q  be  the  number  of  coordinate  pairs,  and 

A;-  be  zero  for  all  line  segments  where  Xx-  >  XSEA  and  x  •+]  <  XLAND. 

The  area  of  the  profile  within  the  provided  boundaries: 

q-1 

A  =  :  A.  (B-1 1 ) 

•  i  £- 


Letting 

C  be  the  appropriate  factor  to  convert  squared  output  units  to  cubed  units. 
The  unit  volume  of  the  profile  within  the  provided  boundaries: 

V  =  AC  (B-13) 


X.  COMPUTING  THE  MEAN  UNIT  VOLUME 

In  the  analysis  module  BEACH,  the  mean  unit  volume  is  computed  for  the  sur¬ 
veys  taken  at  each  profile  line.  If  it  has  been  requested  that  the  unit  volume 
below  vertical  datum  computation  be  performed,  there  is  a  mean  computed  for  the 
unit  volume  below  the  zero  elevation  as  well  as  for  the  part  above.  Otherwise, 
the  mean  unit  volume  is  only  for  that  part  of  the  profile  line  above  the  zero 
elevation.  The  mean  is  a  simple  mean  computed  by  adding  the  appropriate  unit 
volume  during  each  survey  of  the  profile  line  and  dividing  by  the  number  of 
surveys  for  which  the  unit  volume  could  be  computed.  As  discussed  earlier  on 
boundary  determination,  conditions  required  for  allowing  a  computation  of  unit 
volume  below  vertical  datum  are  more  stringent  than  those  required  for  unit 
volume  above.  It  is  thus  possible  that  the  denominator  in  the  computation  of 
the  means  vary.  Using  the  appropriate  number  of  surveys  for  k  and  the  unit 
volume  rather  than  the  distance  xe,  the  computation  for  mean  unit  volume  i;- 
represented  by  equation  (B-4). 

XI.  COMPUTING  STATISTICAL  TRENDS  AND  CORRELATIONS 

Some  basic  statistical  analysis  is  performed  in  the  analysis  modules  BEACH, 
SURVY2,  and  VOLCTR. 
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1 .  Correlation  Coefficients. 


In  analysis  module  SURVY2 ,  the  correlation  coefficient  for  change  in 
shoreline  position  from  the  reference  position  (dependent,  Y)  versus  elapsed 
time  (independent,  X)  is  computed,  if  requested,  for  all  surveys  at  each 
profile  line. 


In  analysis  module  BEACH,  the  following  correlation  coefficients  are 
computed  for  surveys  at  each  profile  line,  if  requested: 

(a)  Change  in  unit  volume  above  the  vertical  datum  ( dependent,  Y) 
versus  change  in  distance  to  shoreline  (independent,  X)  between  consecutive 
surveys . 


(b)  Change  in  unit  volume  above  the  vertical  datum  (dependent,  Y) 
versus  change  in  distance  to  shoreline  (independent,  X)  from  the  selected 
reference  volume  or  position. 

(c)  Change  in  unit  volume  above  the  vertical  datum  (dependent,  Y) 
versus  change  in  unit  volume  below  the  vertical  datum  (independent,  X) 
between  consecutive  surveys. 


(d)  Change  in  unit  volume  above  vertical  datum  ( dependent,  Y )  versus 
change  in  unit  volume  below  the  vertical  datum  (independent,  X)  from  the 
selected  reference  volumes. 


Letting 

k  be  the  number  of  surveys  at  the  profile  line  for  which  both  X  and  Y  can 
be  computed,  the  correlation  coefficient  (Crow,  Davis,  and  Maxfield, 
I9604)  is 


_ kEXY  -  1XEY _ 

>'[kXX2  -  (•  X) "’  ]  I  K  .Y/  -  (EY)‘]  (B-13) 


2.  Regression  Line — Least  Squares  Fit. 

The  regression  coefficient  and  the  regression  line  are  computed  for  the 
following: 

(a)  Change  in  distance  to  shoreline  position  from  the  selected 
reference  position  (dependent,  Y)  versus  elapsed  time  (independent,  X). 

This  computation  is  performed  in  the  analysis  modules  SURVY2  and  BEACH. 

(b)  Change  in  unit  volume  above  datum  from  the  selected  reference 

volume  (dependent,  Y)  versus  elapsed  time  (independent,  X).  This  computa¬ 
tion  is  performed  in  analysis  module  BEACH  only. 


4CROW,  E.L.,  DAVIS,  F.A.,  and  MAXFIELD,  M.W. ,  'tatirtia.n  Manual  ,  Dover 
Publications,  New  York,  I960,  pp.  152  and  158. 
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(c)  Change  in  unit  volume  below  datum  from  the  selected  reference 

volume  (dependent,  Y)  versus  elapsed  time  (independent,  X).  This  computa¬ 
tion  is  performed  in  analysis  module  BEACH  only. 

Letting 

k  be  the  number  of  surveys  of  the  profile  line  for  which  both  X  and  Y 
could  be  computed,  the  regression  coefficient,  slope  of  the  regression 
line  (Crow,  Davis,  and  Maxfield,  I9605),  is 

_  kXXY  -  EXEY  (B-14) 

k£X2  -  (XX)2 

and  the  intercept  (Crow,  Davis,  and  Maxfield,  1960^)  is 


a 


(B— 15) 


3.  Standard  Deviation. 

In  analysis  module  VOLCTR,  the  standard  deviation  for  unit  volume  changes 
within  specific  contour  segments  is  computed.  This  computation  is  used  to 
examine  changes  between  consecutive  surveys  of  profile  lines  at  a  locality  as 

shown  on  TABLE11  (Vol.  I,  Fig.  38)  outputs.  For  profile  lines  surveyed  during 
both  surveys  being,  considered,  the  standard  deviation  is  computed  for  the 
following: 

(a)  Unit  volume  within  each  contour  segment  at  all  the  profile  lines. 

(b)  Total  positive  change  in  unit  volume  at  all  profile  lines.  (Sum 

of  the  positive  changes  in  unit  volume  within  the  contour  segments  at  each 

profile  line  is  the  total  positive  change  at  that  profile  line.) 

(c)  Total  negative  change  in  unit  volume  at  all  profile  lines.  (Sum 

of  the  negative  changes  in  unit  volume  within  the  contour  segments  at  each 

profile  line  is  the  total  negative  change  at  that  profile  line.) 

(d)  Total  change  in  unit  volume  at  all  profile  lines.  (Sum  of  total 
positive  and  negative  changes  at  each  profile  line  is  the  total  change  at 
that  profile  line.) 

Letting 

m  be  a  consecutive  integer  identifying  each  profile  line, 

n  be  the  number  of  profile  lines  at  which  xm  can  be  computed, 


Crow , 

Davis, 

and 

Maxfield, 

op.  cit . ,  p. 

33. 

Crow, 

Davis, 

and 

Maxfield, 

op.  cit.,  p. 

33. 
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X , 


be  the  appropriate  unit  volume  change  at  profile  lines,  and 
be  zero  for  profile  lines  at  which  cannot  be  computed. 


XII.  ELAPSED  TIME 

Elapsed  time,  used  for  both  computational  and  display  purposes  (Vol .  I,  Fig . 
l4),  is  computed  in  analysis  modules  SUPVY2,  PEACH,  and  ELVDIS.  It  may  be  ex¬ 
pressed  as  hours,  days,  months,  or  years  at  the  user's  option,  but  once  the  type 
of  time  has  been  selected,  all  computations  and  displays  for  a  particular  run 
must  be  expressed  in  terms  of  that  time.  No  provision  is  made  for  a  change  ir. 
century;  it  is  assumed  all  data  are  collected  in  the  20th  century. 


Hours. 


or  Years. 


When  elapsed  time  is  expressed  as  hours,  days,  or  years,  the  same 
representation  is  used  both  for  computational  and  display  purposes.  The 
earlier  time  in  the  range  of  dates  covered  by  the  data  is  the  zero  time.  The 
number  of  hours  passed  since  0000  hours  on  1  January  1900  to  the  zero  time  is 
computed.  Next,  the  number  of  hours  passed  since  0000  hours  on  1  January  1900 
to  the  time  of  the  survey  is  computed  and  the  hours  to  zero  time  subtracted 
from  this.  These  hours  become  the  elapsed  time  or  the  hours  are  optionally 
converted  to  days  or  years  if  required. 


The  procedure  is  as  follows: 

(a)  A  12-element  array  containing  number  of  hours  passed  in  a  year 
(excluding  leap  year)  prior  to  the  beginning  of  each  month  in  the  year  is 
established: 


Hx  =  0  H4  =  2,160  H7  =  11,344  Hx  c  =  6,552 

=  744  H5  =  2,880  HP  =  5,088  HX1  =  7,296 

H,  =  1,416  Hp  =  3-624  Ha  =  5,832  Hlr  =  8,016 

(b)  If  the  year  under  consideration  is  a  leap  year  (if  there  is  no 
remainder  when  the  year  is  divided  by  4),  24  hours  are  added  to  each  H 
element  from  on,  inclusive. 

(c)  Unless  the  year  under  consideration  is  1900,  the  number  of  leap 
years  which  have  preceded  the  one  under  consideration  is  computed  by 
subtracting  1  from  the  year  and  dividing  by  4.  The  integer  part  of  the 
quotient  plus  1  is  the  number  of  preceding  leap  years.  If  the  year  is 
1900,  the  number  of  preceding  leap  years  is  zero. 


(d)  Convert  number  of  days  from  preceding  leap  years  to  hours  by 
multiplying  by  24. 

(e)  Find  total  hours  elapsed  by  multiplying  the  year  by  8,760  and 
adding  number  of  preceding  leap  hours.  Add  the  H  element  corresponding  to 
the  month  of  the  survey.  Subtract  1  from  the  day  of  the  survey  and 
multiply  by  24.  Add  the  hour  of  the  survey.  Divide  the  minutes  by  60  and 
add.  Call  this  T  if  the  zero  time  is  being  computed  and  H.  ,f  if  the 
hours  elapsed  to  a  particular  survey  are  being  computed. 


Elapsed  time  in  hours  is  H t  +  -  T, 

Elapsed  time  in  days  is  (H +ct  -  T  )/24 
Elapsed  time  in  years  is  (fi^.  ,,  -  T.)/8,76o 


2.  Months . 


Elapsed  months  is  expressed,  for  display  purposes,  as  the  months  passed 
since  the  beginning  of  each  year.  For  the  graphical  display  of  this  data, 
there  is  a  different  plot  drawn  for  each  year  and  on  each  plot  the  zero  time  is 

1  January  at  0000  of  that  year  (Vol.  I,  Fig.  14).  The  number  of  months  passed 
during  the  year  is  computed  by  subtracting  1  from  the  month  of  the  survey, 
dividing  the  day  of  the  survey  by  30.5  and  adding  the  two  results.  When  the 
day  of  the  survey  exceeds  30  it  is  set  to  30.5  for  the  computation. 

For  computational  purposes,  months  elapsed  is  expressed  as  the  number  of 
months  elapsed  since  the  zero  time: 

(a)  Compute  months  elapsed  during  the  current  year  as  above. 

(b)  Compute  months  elapsed  before  the  current  year  by  subtracting 
the  year  of  the  zero  time  from  the  year  of  the  survey  and  multiplying  by 
12. 


(c)  Add  the  results  of  (a)  and  (b)  to  obtain  total  elapsed  months. 
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APPENDIX  C 


BP AS  ERROR  MESSAGES 

The  errors  in  the  BPAS  are  in  two  categories — those  which  are  informative 

and  those  which  are  fatal.  Each  error  has  been  assigned  a  unique  number  which 
is  printed  along  with  the  error  messages. 

The  informative  messages,  identified  by  an  "I,"  do  not  always  flag  a  real 
error  and  have  been  included  to  alert  the  user  that  the  program  has  encountered 
an  unusual  or  potentially  serious  situation. 

The  fatal  errors,  identified  with  an  nF,"  cause  the  programs  to  stop 
execution  after  they  are  printed.  These  errors  occur  when  further  execution  is 
impossible  or  it  is  probable  that  further  execution  will  not  produce  usable 
results. 

Table  C-l,  which  is  a  sequential  list  of  the  BPAS  error  messages  by  number, 
contains  the  error  message  which  is  printed,  the  category  of  the  error,  and  the 
action  that  should  be  taken.  The  remarks  column  contains  further  information, 
such  as  the  most  common  sources  of  the  problem  and  steps  which  may  help  circum¬ 
vent  the  problem.  Some  errors  can  be  eliminated  only  with  a  program  modifica¬ 
tion;  however,  a  discussion  describing  how  to  implement  such  modifications  is 
beyond  the  scope  of  this  report. 
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•“.able  C-l.  Listing  of  BPAS  error  messages. — Continued 

Error  Ho.  Description  Severity  Action  Remarks 


O  S-  C  T? 
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able  0-1  .  l.ijtiru'  of  i'PAC  error  metis ares. — 0^ntinue<i 
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These  errors  suggest  processing  through  analysis  routines  of  data  which  have  not  been  created  directly  from  the  EDIT? 


Table  C-l.  Listing  of  BPAS  error  messages. — Continued 
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No  consecutive  surveys  of  line  —  I  None  required.  See  comments  Tor  error  WJ. 

had  a  datum  intercept .  No  output 
for  this  line. 


Table  C-l.  Listing  of  BPAS  error  messages. — Continued 

top  No.  Description  Severity  Action  Remarks 
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These  errors  suggest  processing  through  analysis  routines  of  data  which  have  not  been  created  directly  from  the  ED1T2 
program. 
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(position  69)  must  be  1  or  0.  F  Correct  and  PLOT  spec  card  improperly  completed, 

red  for  plot  type - -  retry. 


Table  C-l.  Listing  of  BPAS  error  messages. — Continued 

Error  No.  Description  Severity  Action  Remarks 


i 


defined  for  statistical  correlations.  retry  if  blank 

TAPLE8A  suppressed.  TABLE8A  output 

is  desired. 


APPENDIX  D 


DATA  REQUIREMENTS 


1 .  Basic  Assumptions. 

The  data  for  which  the  BPAS  was  designed  consist  of  beach  profile  data  and 
the  computations  performed  deal  mainly  with  changes  in  shoreline  and  other 
contour  positions  and  unit  volumes.  Figure  D-1  is  a  schematic  of  a  typical 
beach  profile  survey.  The  following  restrictions  are  placed  on  the  data  which 
can  be  processed  by  the  system: 

(a)  Distances  to  each  successive  surveyed  point  must  be  greater  than 
or  equal  to  the  distance  to  the  previous  one. 

(b)  Each  survey  must  begin  landward  of  the  shoreline  position  and 
proceed  seaward. 

(c)  There  must  be  at  least  one  point  with  an  elevation  greater  than 

zero. 

Based  on  these  restrictions  there  are  three  types  of  profiles  which  can  be 
analyzed  by  the  system  (see  Vol.  I,  Sec.  IV ): 

(a)  Type  1:  The  profile  extends  from  a  point  landward  to  a  point 
seaward  of  the  shoreline. 

(b)  Type  2:  The  profile  extends  only  to  the  shoreline  or  the 
shoreline  position  can  be  extrapolated. 

(c)  Type  3-  The  profile  does  not  extend  to  the  shoreline  and  the 
shoreline  position  cannot  be  extrapolated. 
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2.  Program  Limitations. 


Because  of  array  dimension  sizes  and  other  computational  procedures,  the 
following  additional  restrictions  are  placed  on  the  data: 

(a)  There  must  be  at  least  two  but  no  more  than  60  surveyed  points 
defining  a  profile. 

(b)  Profile  line  and  survey  numbers  must  be  numeric.  The  difference 
between  the  largest  and  smallest  profile  number  must  not  exceed  99;  between 
the  survey  numbers,  no  more  than  149.  Profile  line  numbers  must  fall 
between  1  and  999,  survey  numbers  between  1  and  9999. 

(c)  No  more  than  50  years  may  be  represented  in  a  single  data  set. 

(d)  All  data  on  a  final  data  set,  ready  for  input  to  the  analysis 
routines,  should  be  in  the  same  units  and  referenced  to  the  same  vertical 

and  horizontal  datum. 

(e)  In  computations  using  the  distance  between  profile  lines,  the 
lines  are  assumed  to  be  parallel. 
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APPENDIX  E 
SAMPLE  CODING  FORMS 
FOR 

SURVEY  INPUT  DATA 
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PROFILE  SURVEY  DATA 

CODING  FORM 


BEACH  EVALUATION  PROGRAM 

SURVEY  REPORT  coastal  engineering  research  center 

(  INITIAL  SHEET  )  KINGMAN  BUILDING 

FT.  BELVOIR,  VA  22060 


LOCALITY  CODE 


NAME  OF  BEACH 


12  3  4  5 


SURVEY  NUMBER  111  I  I 

6  7  8  9 

SHEET  NUMBER  □□  OF  □□ 

10  It  i 2  i 3 

v  R  MQ  DAY 

SURVEY  DATE  QH  Qj  [JJ 

14  15  16  17  18  19 

TIME  CD  m 

22  23  24  25 


DATUM  □ 

26 


0  SWL  3  MTL  6  MHW  9  MLL 

1  NGVD  4  MLW  7  IGLO 

2  MSL  5  MLLW  8  LWD 

A  OTHER  _ 


DATA  IS  COLLECTED  IN  fTI 

27  28 

(  FT,  M,  ETC  ) 


LOCATION  l  PROFILE  OR  PROFILE  NO.) 

SURVEY  METHOD 

1  LEVEl  AND  TAPE 

2  STADIA 

3  EMERY  Q 

4  HAND  LEVELING 

5  OTHER 

IF  THERE  HAS  BEEN  ANY  CHANGE 
IN  PROFILE  BENCHMARK  SINCE 
LAST  REPORT, CHECK  BOX  (/) 

AND  COMMENT  ON  BACK  OF  THIS 
FORM  □ 

IF  YOU  HAVE  ANY  OTHER  COMMENT 
CHECK  BOX  ;/)  AND  USE  THE  BACK 
OF  THIS  FORM  Q 


DISTANCE  MUST  BE 

IN  ASCENDING  ORDER  ELEVATION 


4- 

□ 

C  1000)0001  f  1 01  (It 

rrm 

4- 

□ 

(ioohioi  m 

rm 

I.H 

□ 

31 

♦ 

32  33  34  35 

36 

4- 

37  38  39 

40 

□ 

Mill 

n 

rm 

□ 

4  1 
+ 

42  43  44  45 

46 

4- 

47  48  49 

50 

□ 

1  1  1 

□ 

rm 

□ 

5i 

+ 

52  53  54  55 

56 

57  58  59 

60 

□ 

□  1  1  1 

Q 

rm 

□ 

61 

■f 

□ 

62  63  64  65 

66 

4- 

67  68  69 

70 

mil 

□ 

rm 

□ 

7  1 

n  73  74  75 

76 

7/  '8  79 

80 

PARTY  CHIEF 


FORM  FILLED  87  DATE 


profile  survey  data 

CODING  FORM 


SURVEY  REPORT 

( CONTINUATION  SHEET) 


BEACH  EVALUATION  PROGRAM 
COASTAL  ENGINEERING  RESEARCH  CENTER 
KINGMAN  BUILDING 
FT.  BELVOIR,  VA.  22060 


LOCALITY  CODE 


12  3  4  5 


SURVEY  NUMBER  LI 


6  7  8  9 


SHEET  NUMBER  QH  OF  LD 

10  11  1 2  i 3 


NAME  OF  BEACH 


location  (profile  or  profile  no.) 


DISTANCE 

(MUST  BE  IN  ASCENDING  ORDER) 


ELEVATION 


IIOOOIOOOHIOI  (II 

•f 

(1001(10)  III 

Ul 

□ 

rrm 

□ 

rr 

□ 

□ 

21 

22  23  24  25 

26 

27  28 

29 

30 

+ 

+ 

□ 

□ 

LQ.-1J 

□ 

CE 

□ 

3i 

32  33  34  35 

36 

37  38 

39 

40 

+ 

4- 

□ 

□ 

rrm 

□ 

LL 

□ 

41 

42  43  44  45 

46 

47  48  49 

50 

+ 

+ 

□ 

□ 

rrm 

□ 

LU 

□ 

5) 

52  53  54  55 

56 

57  58 

59 

60 

■f 

4- 

□ 

□ 

Mill 

□ 

or 

□ 

si 

62  63  64  65 

66 

67  68 

69 

70 

4- 

4- 

□ 

□ 

rrm 

□ 

m 

□ 

7 1 

72  73  74  75 

76 

77  78 

79 

80 

ATTEND  IX  F 


FAKPLE  SURVEY  IIirUT  DATA,  EDIT1  AND  EDITS 

These  sample  input  data,  input  to  the  EDIT1  and  EDITS  programs; ,  were  used 
to  create  the  final  data  file  used  by  the  analysis  routines.  All  sample  cut- 
puts  in  the  series,  except  those  which  demonstrate  features  net  availalle  from 
these  data,  were  created  using  the  input.  Any  output  not  created  by  these  da* 
is  so  identified. 
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APPENDIX  C. 


E AMPLE  FINAL  DATA  FILE,  EDIT] 

These  data  were  created  by  the  FL’ITL  program  fr'  n.  the  samp]  e  input  data  in 
Appendix  F.  Leeause  data  at  profile  lines  68  ar.d  69  were  collected  tc  a 
different  vertical  datum  (pier  deck)  than  at  other  profile  lines,  a  correction 
of  feet  was  made  t o  each  vertical  cc ordinate  at  lines  66  and  69.  These 

data  are  in  the  card  irr.are  format  to  make  them  more  legible.  These  data,  in 
magnetic  media  format,  were  used  to  produce  all  the  analysis  module  sample 
outputs  used  this  series  of  reports,  except  those  otherwise  identified. 
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